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In terms of Bogoliubov-de Gennes approach, we investigate Majorana bound state (MBS) in vortex 
of proximity-induced superconductivity on the surface of topological insulator (TI). Mapping out the 
local density of states (LDOS) of quasiparticle excitations as a function of energy and distance from 
vortex center, it is found that the spectral distribution evolves from “V”-shape to “Y”-shape with 
emergence of MBS upon variation of chemical potential, consistent with the STM/STS measurement 
in a very recent experiment [Xu et al, Phys. Rev. Lett. 114, 017001 (2015)] on Bi 2 Te 3 thin layer on 
the top of NbSe 2 . Moreover, we demonstrate that there is a “checkerboard” pattern in the relative 
LDOS between spin up and down channels, which maps out directly the quantum mechanical wave 
function of MBS. Therefore, spin-resolved STM/STS technique is expected to be able to provide 
phase sensitive evidence for MBS in vortex core of topological superconductor. 

PACS numbers: 03.65.Vf 74.25.Ha 74.45.+C 74.55.+V 03.67.Lx 


Introduction .— Majorana bound states (MBSs) are 
under intensive search in topolgoical superconductors 
(SCs) H i. Because of the peculiar property that a 
particle is equivalent to its antiparticle [3|, non-Abelian 
quantum statistics can be generated which is believed to 
be potentially important for achievin g de coherence-free 
topological quantum computation |4l4lQj. MBSs were 
first predicted in the spinless pairing states of fermions 
including the Pfaffian state of quantum Hall system 
with 5/2 filling [ll| and ultracold atomic gases with p- 
wave Feshbach resonance [l2[ . Recently, heterostructures 
made of 5-wave SC and topological insulator (TI) [isl, or 
the combination of semiconductor with Rashba spin-orbit 
coupling (SOC) and ferromagnetic insulator |14l4l6|. are 
proposed for realizing topological SC with the spin de gree 
of freedom suppressed by SOC. Sr 2 Ru 04 [I^, UPtq ]l8| 
and doped TI Cua,Bi 2 Se 3 0 are also discussed as pos¬ 
sible candidates. 

MBS is a unique quasiparticle excitation of topological 
SC carrying zero energy and zero total angular momen¬ 
tum, which appears at ends of one-dimensional systems 
or at vortex cores, where the SC gap is closed |13I416U 2Q|- 
[Hl . The unique property of equivalence between particle 
and antiparticle, which renders MBS noble for advanced 
app lications, makes them difficult to be identified HHi- 
Is if . The hurdle has been tackled experimentally. Up to 
this moment, there are already several reports on possi¬ 
ble MBSs in ID nanowires 32-[35j. where quasiparticle 
states are observed at ends of nanowires at zero-energy 
bias. 


Very recently scanning tunneling microscopy and spec¬ 
troscopy (STM/STS) experiments have been performed 
in vortex state of proximity-induced SC in Bi 2 Te 3 thin 
layer [36|. Mapping the local density of states (LDOS) 
measured by differential conductance dljdV as a func¬ 
tion of bias voltage and the distance from vortex core. 



FIG. 1: Schematics for the model system in the present 

study, (a) System geometry with a topological insulator (TI) 
slab on top of superconductor (SG) carrying on one vortex, 
with the red dots for Majorana bound states, (b) Dispersion 
relation of TI with helical surface states (in blue color) and 
the bulk valence and conduction bands (in gold color) with 
a gap of 2eo. The chemical potential p is measured from the 
Dirac point of the topological surface states. 


an evolution from “V”-shape for thin TI to “Y”-shape 
for thick TI was observed, and it is conjectured that the 
“Y”-shaped LDOS with its center at zero-energy bias is 
due to the appearance of MBSs. 

In this work, we examine the same geometry as in the 
experiment [36| (see Fig. [T]) in terms of Bogoliubov-de 
Gennes (BdG) analysis based on a model system includ¬ 
ing explicitly the topological surface states of TI slab with 
proximity-induced SC. We clarify that, due to the rela¬ 
tion among energy, angular momentum and spatial dis¬ 
tribution of quasiparticle excitation, the spatial-energy 
distribution of LDOS evolves from “V”-shape to “Y”- 
shape corresponding to the absence and presence of MBS 
at vortex core when the chemical potential is reduced, 
which is in nice agreement with the STM/STS measure¬ 
ment [ 3 ^. Moreover, exploring the total angular mo¬ 
mentum of quasiparticle excitation contributed from spin 
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and orbital angular momenta and the phase winding of 
vortex, we demonstrate that there should be a “checker¬ 
board” pattern in spin-resolved LDOS associated with 
MBS. This suggests a phase sensitive way for identifying 
MBS as a single quasiparticle by spin-resolved STM/STS 
technique. 


Model Hamiltonian and BdG approach .— As can be 
read from Fig. 4 of Ref. [s^, when the thickness of TI 
slab increases, the Fermi level shifts from a position cut¬ 
ting the bulk conduction band [for one to four quintuple 
layers (QLs)] to a position below the conduction band 
(for five and six QLs), while the size of band gap be¬ 
tween the bulk conduction and valance bands remains 
almost unchanged. For the whole range of thickness, the 
SC gap can be taken as constant as a fairly good ap¬ 
proximation (see Fig. 2 of Ref. Q). In order to model 
the situation realized in the experiment, we consider a 
system of TI slab carrying on topological surface states 
and with proximity-induced SC as shown schematically 
in Fig. [IJa), where the chemical potential is varied while 
the thickness of TI slab and the SC gap are fixed. The 
BdG Hamiltonian is given by 


^BdG = 

with the effective Hamiltonian of TI 


HTiir) A(r) 
At(r) 


Hti — — ivF 

and Dirac mass 
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in the cylindric coordinates, where Si and di are Pauli 
matrices for spin up and down states and orbitals of Bi 
and Te with even and odd parities respectively. The bulk 
gap between valance and conduction bands is taken as 
2eo and Fermi level /i is measured from the Dirac point 
of the topological surface states of TI; /i > eo corresponds 
to thin TI slabs and < eo to thick TI slabs in the exper¬ 
iment [ 3 ^. Since the energy gap of Bi 2 Te 3 is eo — O.leV 
and the proximity-induced SC gap is Ao — 1 meV in 
the experiment [36|, we take Ao/eo = 0.02 for numer¬ 
ical calculation in the present work, noticing that this 
value influences substantially the distribution of spectral 
weight. The Fermi velocity of topological surface state 
corresponds the slope of linear dispersion in Fig. [11(b), 
which yields the unit for lateral length = v-f/cq. 
Taking vf — 0.2 nm-eV for Bi 2 Te 3 from another exper¬ 
iment one has 2nm. The effective mass is 

taken as m = 2eo/v^ for simplicity. 

For simplicity we consider the case with only one vor¬ 
tex, where the SC gap takes the form A = e'^^AoiSy with 
(j) the azimuthal angle in the cylindrical coordinate. The 



FIG. 2: (a) and (c) Eigen energy of quasiparticle excitation 
in vortex as a function of total angular moment j. The red 
dot in (a) indicates two degenerate MBSs localized at the top 
and bottom surface of TI slab. Colors in (c) are for values of 
momentum along z direction, (b) and (d) LDOS in a window 
close to zero energy |F|/eo < 0.002 as denoted by the gray bar 
in (a) and (c). Parameters are Ao/eo = 0.02 and t = Gk^^. 


eigen wave functions of quasiparticle excitations take the 
form 


f 

V 


gi-U l + D’i’Uj ^{r, z)'\ 


where Uj e ^^e 5 electron and hole 

wave functions, are two-component vectors referring to 
the two orbitals, and j and E are the angular momen¬ 
tum and eigen energy. For numerical diagonalization, 
we expand the wave functions Ujj^{r) and in 

terms of the Bessel orthonormal functions in radial di¬ 
rection [sil and the Gauss-Lobatto functions in z direc¬ 
tion [ 3 ^. Since the coherence length ^ = vf/Aq = bOk^^ 
is large, a radius R = 2b0kQ^ is used for numerical cal¬ 
culations. 

When the Fermi level lies in between the bulk conduc¬ 
tion and valence bands, namely ja < cq (thick TI slab 
in experiments), SG gap is opened only in the surface 
Dirac dispersions of TI. As the result, two MBSs with 
zero energy and zero angular momentum appear at the 
top and bottom surfaces of TI slab as shown in Fig. O 
The energy difference between MBSs and the lowest ex¬ 
citations is 6E 0.2Ao(^ SAg/yu) for yu/e = 0.5. The 
spatial distributions of the two MBSs on the top and bot¬ 
tom surfaces are the same since in the present simplified 
model the effect of SG substrate is taken into account 
by a uniform proximity-induced SG gap. Because the 
two MBSs are localized at the two surfaces of TI slab. 
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FIG. 3: (a)-(d) LDOS as a function of radius r and energy E 
with smearing factor ?7/eo = 0.004. The spectra are taken at 
z = 0.24/c^^ from the surface of the TI slab. Dotted curves 
are for eye guide, (e) Spatial distribution of MBS, with open 
circles for numerical results and solid curves for analytical 
results \uq^q\ = Jo{rkF) and |'Uo,ol “ Ji{rkF) 
where ^ = Aq/vf = 50/c^^. All parameters are the same 
as in Fig. [2 


they are well described by the 2D model 


analytical solution [25 


and the 


When the Fermi level falls into the bulk conduction 
band, namely /i > cq (thin TI slab in experiments), low- 
energy quasiparticle excitations carrying finite momenta 
kz = nir/t with t the thickness of TI slab spread along 
2 ; direction in vortex core as can be seen in Figs. He) 
and (d). Due to these states in the normal vortex core 
the two MBSs on surfaces interact with each other [ 2 ^ 
yielding a mini gap as can be seen in inset of Fig. He). 

LDOS for quasiparticle excitations are given by the 
electron wave functions. 


which is measured directly by the differential conduc¬ 
tance dl/dV in STM/STS experiments. Judging from 
the continuous spectra observed in Fig. 2 in Ref. [s^, 
the energy resolution in experiments is larger than the 
energy differences between quasiparticle levels. In order 
to simulate this situation, we replace 6{x) with a smear¬ 
ing function C{x,r]) = x/[7r{r]‘^ + x‘^)] with r]leo = 0.004 
in Eq. ([5]). 

The spatial distribution of a quasiparticle excitation 
is mainly governed by the orbital angular momentum, 
which appears in the wave functions approximately given 
by Bessel functions. The energy dispersion with respect 
to the total angular momentum then gives a correlation 
between energy and spatial position of maximal abso¬ 
lute value of quasiparticle wave function. As a general 
rule, quasiparticles with higher energy locate at a posi¬ 
tion farer away from the center of vortex. As shown in 
Fig. [3l for /i/eo = 0.5 and 1 (thick TI slab in experi¬ 
ments), accompanied by MBSs with zero energy locat¬ 
ing around the center of vortex core [see Figs. [Sfa) and 
(b)], the LDOS takes a “Y”-shape; for /i/eo = 2 (thin 
TI slab in experiments), the LDOS takes a “V”-shape 
corresponding to the absence of spectral weight at zero 
energy [see Figs. He) and (d)]. This tendency can also be 
seen directly from the wave function of MBSs. As shown 
in Fig. He), the MBSs distribute over a larger area in 
the vortex core when the chemical potential is lowered, 
since their wave function is given by the Bessel function 
Jo{rk-F) [111, [ 2 ^, where /cp, the intercept of Dirac disper¬ 
sion with Fermi level, decreases with decreasing /i upon 
lowering the Fermi level [see Fig. [II[b)]. The above evo¬ 
lution of spectral shape of LDOS is consistent with the 
tendency observed in the STM/STS experiments where 
the thickness of TI slab is changed systematically [s^, 
and supports the existence of MBSs in thick TI slabs. 

LDOS presented in a previous work exhibits asymme- 
try in spectrum with respect to the zero energy [23| . The 
apparent difference between them and the ones in Fig. [3] 
is due to that a large SC gap Aq 10 meV was taken in 
that work as compared with a small smearing factor. 

Spin-resolved LDOS .—MBSs as quasiparticle excita¬ 
tion with zero energy and zero angular momentum can 
be better identified by spin-resolved LDOS even with the 
same STM/STS resolution. To demonstrate this idea, 
we check the wave functions of quasiparticle excitations 
paying attention to spin degree of freedom. Due to the 
rotational symmetry, the total angular momentum j 

j = l + s/2- 1/2 (6) 


Af{E,r,z)=M^iE,r,z)+Af^iE,r,z) (4) 

with 

KiE, r,z) = Y^ \ulE'{r)\^KE - E'), (5) 

E'd 


contributed from orbital and spin angular momenta 
and the phase winding of SC gap (presumed as anti¬ 
clockwise) is conserved. Therefore, it is clear that the 
MBS with zero total angular momentum has two compo¬ 
nents with / = 0 for up spin (5 = 1 ) and / = 1 for down 
spin (5 = —1). The detailed wave functions of MBS are 
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FIG. 4: Wave functions of several low-energy quasiparti¬ 
cle excitations in vortex of topological superconductor. Pa¬ 
rameters are the same as Fig. [2] except for z—^.2AkQ^ and 
/i/eo = 0.5. 


displayed in Fig. [4| with spin up and down separately, 
where the amplitudes oscillate with distance from vortex 
center in terms of Bessel functions 25, [ 2 ^. According 
to our numerical calculations, the excited state with en¬ 
ergy E = O. 2 A 0 and total angular momentum j = —1 
exhibits a spatial distribution almost the same as the 
MBS [see Fig. [4]. It is intriguing to notice that, how¬ 
ever, the oscillations in the two wave functions are out of 
phase when the spin state is specified. The same out-of- 
phase oscillations happen in the excited state with energy 

= — O. 2 A 0 and total angular momentum j = 1. 

This property can be used for distinguishing MBSs 
from the first excited states. In order to see this ex¬ 
plicitly, we calculate the spin-resolved LDOSs and 
A/ 4 , given in Eq. ([5]) and evaluate the relative LDOS by 
taking the ratio between them. The result thus obtained 
for /i/eo = 0.5 and Ao/eo = 0.02 is displayed in Fig. [3 
where a checkerboard pattern is observed in the spatial- 
energy mapping of LDOS with the period of 10 nm in 
space and 0.4 meV in energy. We emphasize that, while 
the smearing factor in Figs. [3] and [5] is the same, corre¬ 
sponding to the same STM/STS resolution, the relative 
spin-specific LDOS resolves the wave function of MBS up 
to the quantum limit. In a sharp contrast, the spectrum 
of quasiparticle excitations for ja/eo = 2 is continuous, 
since bound states with kz ^ 0 spread out in z direc¬ 
tion in the vortex core, and the relative spin-resolved 
LDOS loses the oscillating phase. In an anti¬ 

vortex with clockwise phase winding, the last term in 
Eq. ([ 6 ]) should change sign, and it is easy to see that the 
checkerboard pattern in Eig. [3 should switch the purple- 
blue contrast. Therefore, the checkerboard pattern of 
spin-resolved LDOS for quasiparticle excitation displayed 
in Eig. [3 provides a phase sensitive evidence for MBS in 
vortex core of topological SC. It is worth noticing that, if 
the energy resolution of STM/STS measurement is high 
enough, checkerboard patterns can be detected even in 



FIG. 5: Relative spin-resolved LDOS A/f/A/*^ as function of 
energy and distance from vortex center with the same smear¬ 
ing parameter as in Fig. [3 Parameters are the same as Fig. (4] 


individual spin-resolved LDOS (without taking the ratio 
between signals in two spin channels, see Appendix). Ex¬ 
perimentally, spin-resolved LDOS can be obtained by us¬ 
ing a spin-polarized STM tip [i^. It is noticed that, with 
very high resolution splitting in spectral weight due to 
spin states can be observed even in the total LDOS |41|. 


Conclusions .— We clarify that, in a vortex core of 
proximity-induced SC on the surface of TI, the spatial- 
energy mapping of local density of states of quasiparticle 
excitations evolves from “V”- to “Y”-shape with emer¬ 
gence of Majorana bound state, which is in good agree¬ 
ment with a very recent STM/STS experiment. More¬ 
over, we demonstrate that the relative density of states 
between spin up and down channels exhibits a checker¬ 
board pattern reflecting the out of phase oscillations in 
the spin-specific wave functions of Mojarana bound state 
and excited states. It provides a clue toward observing 
Majorana bound state in terms of phase sensitive signal 
by using spin-resolved STM/STS experiments. 
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Culture, Sports, Science and Technology of Japan. 
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Appendix: LDOS in individual spin channels 

Spin-resolved local densities of states (LDOS) shown 
in Fig. [5] are re-plotted separately in two spin channels 
in Fig. [6l Checkerboard observed in the relative spin- 
resolved LDOS cannot be seen clearly in the two individ¬ 
ual spin-resolved LDOS, which demonstrates the advan¬ 
tage of taking their ratio in highlighting the oscillating 
behavior in LDOS associated with the Majorana bound 
state. The LDOS is symmetric (asymmetric) in the spin- 
up (-down) channel corresponding to the anti-clockwise 
phase winding of vortex, as can be understood from the 
wave functions exhibited in Fig. [H 

Meanwhile, it is worth noticing that checkerboard pat¬ 
terns can be identified in LDOS in individual spin chan¬ 
nels provided the energy resolution of STM/STS mea¬ 
surement is high enough. In order to demonstrate this 
point, we simulate a case with smearing factor = 
0.002 smaller than the one for Fig. [6] ( 77 /eo = 0.004) and 
display the results in Fig. 0 Here, once again checker¬ 
board pattern is enhanced in relative LDOS. 



f;/Ao e/Ao f;/Ao 

FIG. 6 : LDOS in the two spin channels as a function of ra¬ 
dius r and energy E. Parameters are the same as in Fig. [S] 
For comparison, the relative spin-resolved LDOS is displayed 
which is the central part of that in Fig. O 
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